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Abstract: If dark matter (DM) annihilation accounts for the tantahzing excess of cosmic 
ray electron/positrons, as reported by the PAMELA, ATIC, HESS and FERMI observa- 
tories, then the imphed annihilation cross section must be relatively large. This results, in 
the context of standard cosmological models, in very small relic DM abundances that are 
incompatible with astrophysical observations. We explore possible resolutions to this ap- 
parent conflict in terms of non-standard cosmological scenarios; plausibly allowing for large 
cross sections, while maintaining relic abundances in accord with current observations. 



Keywords: Low scale B - L, Dark Matter 



Contents 

|l]. Introduction ^ 

|2|. PAMELA, ATIC, FERMI and HESS Anomalies | 

|3|. Relic Abundance and Non-conventional Cosmology ^ 

^ Extra Constraints on the Annihilation Cross Section and the Viability 

of our Model ^ 

[H]. Conclusions 10 



1. Introduction 

The context of our investigation is set by the overwhelming evidence for a non-baryonic 
component dominating the matter content of our Universe, as inferred from a combination 
of colluding astrophysical phenomena (the dynamics of galaxies and clusters; large scale 
structure in the Universe; cosmic microwave background fluctuations; big bang nucleosyn- 
thesis). Though relying exclusively on large scale gravitational signatures, the implied 
existence of DM alludes to microscopic fundamental physics beyond the standard model. 
But although there are several ongoing attempts at direct detection of expected DM candi- 
dates, apart from the controversial results of DAMA collaboration [Q] and recent tentative 
data points emanating from the CDMSII experiment their outcomes have hitherto 
been invariably negative. It is in this context that the significant excitement concerning 
the recent measurements of excess electron and positron flux, as reported by the PAMELA, 
ATIC, HESS and FERMI collaborations, among others, arose. 

PAMELA'S observations, reporting excess flux between 8 and 80 GeV Q, with no 
excess in the corresponding anti-proton flux, confirm and extend previous results obtained 
by HEAT g] and AIMS g. The ATIC @ experiment data, on the other hand, show 
significant excess electron and positron flux at energies around 300 — 800 GeV. 

Although the PAMELA data is marginally consistent with calculations of the cosmic 
ray background employing a 'soft' background electron fiux spectrum the ATIC and 
FERMI observations more clearly point to the existence of an additional source of high 
energy electrons and positrons; and since such particles cannot travel very far without 
much energy loss, their source must be local ( < 1 kpc from the solar system). There 
are plausible astrophysical explanations for this excess too, e.g. in terms of local pulsars 
and supernovas remnants but they seem less natural than when invoked in the low 

energy region of the spectrum. The excesses in flux could also result from DM annihilation 
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or decay, and a vast literature has recently arisen around the subject (see, e.g., 10] for 
a review), partly instigated by a salient feature (perhaps especially prominent in the ATIC 
measurements) of the cosmic ray positron spectrum produced by DM annihilation; the fact 
that it drops off at Ef,+ = m^, while the flux produced by, e.g., a single pulsars falls off more 
gradually. Though this feature seems significantly more subdued in the FERMI data |2|, 



the relevant region is relatively less well represented there — at least in comparison with 
the HESS data, where a palpable drop-off is present |13|, [14|. Nevertheless, is difficult 



to conclusively differentiate between astrophysical and DM descriptions of the observed 
positron excesses in any conclusive manner, although future experiments may be able to 
achieve this. One advantage of DM models is that, despite their diversity, their cosmic 
rays emanate from far simpler physical objects and are therefore far easier to falsify than 
those pertaining to astrophysical sources. 

In any case, if the DM scenario is to explain the observed anomalous flux, one major 
hurdle has to be surpassed. At present the major difficulty facing the DM annihilation 
is the large cross section apparently required to fit the excess flux — one that seems 
incompatible with straightforward estimates of the relic DM abundance in conventional 
cosmological models. 

This paper is an attempt at reconciling the aforementioned observed cosmic ray excess 
with an interpretation of the results in terms of annihilating halo dark matter particles. 
In the next section we argue that a widely disseminated solutions to this problem, in 
terms of boost factors, including those sprouting from Sommerfeld enhancement and similar 
effects, seems untenable. Prom this stems our motivation for investigating non-standard 
cosmological models, allowing for larger DM annihilation cross sections while being less 
severely constricted by current observations. 

2. PAMELA, ATIC, FERMI and HESS Anomalies 

The PAMELA experiment measured an excess of cosmic ray positrons with no indication of 
any excess of anti-proton flux. Therefore, if it is indeed the dark matter that is responsible 
for the positron excess, it seems natural to consider a type of DM particle that annihilates 
predominantly into channels. The positron flux in the galactic halo is given in terms 
of the production rate of the positrons from DM annihilation, which is given by 



where {av) f = aj refers to the averaged annihilation cross section into the final state / 
and {dN/dE)f is the fragmentation function, representing the number of positrons with 
energy E, produced from the final state /. p{y) is a DM halo mass profile. Though there 
are several types of proposed halo dark matter density profiles, we adopt here the standard 
NPW profile pl. 



Pig. H shows the total absolute flux, $*<J<^' + $*ot«' ^ generated by the DM annihilation 
into e^e~ , and t~^t~ as function of positron energy for niy = 1 TeV and 10~^ GeV~^. 
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Figure 1: The total absolute flux in units of GeV^m^^s~-'^sr~^, generated by the DM annihilation 
into e^e^ (green line), /i'^/x" (black line) and r"'"T~ (red line), as function of positron energy for 
— 1 TeV, for thermal averaging cross section fO~^ GeV~^. 



In our analysis, the MED diffusion model of Delahaye et. al. [16] and NFW|I5|] Galactic 



halo with scale-length 20 kpc are assumed. Also, <|)*'^"' and <I)*°*'^' are defined as 

f^total ^ ^mi _^ ^sec ^ ^total ^ ^prim _^ ^sec ^2.2) 



The flux of positrons<I>^ is given from the number density of positron through several 
steps(for instance, see ref. [^). 

The fragmentation function for direct process is almost monotonic; (4^) ~ S{E — 



' ee 

c)- For {%)^^ and (^)^^, we referred to the ref. @. 



Calculated in terms of fitting functions matching the fluxes deduced via standard sim- 
ulations of cosmic ray production and propagation [20|, the astrophysical background 



fluxes of positrons; ^^^"^ and $^1'^, are given by the followings: 

$^"™-(e) = ,,J cm-h-hr-^), 

^ ^ 1 + IleO-9 + 3.2e2-i5^ " 

'-^'> 1 + 11061-5 + 600e2-9 + 580e4-2^'^"' ^' 

Kfi^) = j^(cm"2s-isr-i), (2.3) 

^ ' 1 + 650e2-3 + l500e4-2 ^ ^' ^ ^ 

where e = E/{lGeV). 

We use a DM particle of mass 1 TeV, consistent with the maximum allowed by FERMI 
7-ray observations of dark matter dominated virilized gravitational systems |^T|, [^2|. As 
is immediately apparent, and are quite small at all energies, and can be neglected 
with respect to ^p^*"^; the same goes for the computed positron flux for (av) ~ 10~^ GeV~^, 
which remains far below the background. Therefore, in such a situation, and in terms of 
the plotted quantities, any excess flux would have to be explained in terms of modifications 
to the astrophysical background (including the possible effects of local pulsars, supernovae 
remnants etc.). Alternatively, a larger cross section can lift the computed DM flux above the 
background, so as to explain any apparent excess. It turns out that a (av) ~ 10~^ GeV~^ is 
compatible with PAMELA observations (without invoking any changes in the background 
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model). In addition, for positron energy E > 200 GeV, the positron flux exceeds the 
electron background, which fares well in explaining the ATIC, HESS and FERMI excess 
at 300 - 800 GeV. 

Positron Fraction : Ft; IVI^„=1000 (GeV) and (ov) =10"" (GeV)'^ 

10° ^ ■ ' — ' — ' — ' ■ — ' — ' — ' ■ — ' — ■ — ' ' 




Figure 2: The positron fraction for DM annihilation into e e+ (green line), fi^^ (black line) and 
T+T^ (red line) for dark matter mass = f TeV and thermal averaging cross section fO~^ GeV~^. 

In Fig. we plot the positron fraction 

q>DM ^ ^lf{E) + $^"'"(^;) + ^1<L%E) ' 

the quantity that PAMELA actually measured, for DM annihilation into e~^e^ , ^~ and 
r+r~, along with the relevant measurements. It is clear, as expected, that DM annihila- 
tion with {av) ~ 10-6 GeV-2 can easily account for the PAMELA measurements while 
employing a standard astrophysical background. These measurements explore the exis- 
tence of positron excess for energies confined between 8 and 80 GeV. In the next PAMELA 
experiment, the search for excess flux will be extended up to energies ~ 270 GeV. We note 
here that, at such energies, the positron flux due to DM annihilation is of order the electron 
primary flux and so the expected positron fraction would tend toward 

MI 

: 2l -0(0.3). (2.5) 

\ e~ e~ ' e+ 

Such a result would represent a significant signature confirming that DM annihilation 
is the source of this observed positron excess, making the case for an increased DM flux 
magnitude (compared to what can be inferred from calculations using standard smooth halo 
density profiles and cross sections compatible with thermal relic abundance calculations), 
particularly more pressing. 

Already, myriads of remedies, invoking various 'boost factors', have been proposed. 
Including, for example, in terms of increased DM density inside the galaxy's the subhalo 
population. Yet very little such substructure is expected in the solar neighborhood; efficient 
stripping dissolves subhalos in these inner regions, ensuring that most mass in subhalos is 
concentrated at far larger galactocentric radii. And positrons with energies of order 100 
GeV or above are expected to emanate from DM annihilations within a few hundred pc of 
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the solar system. This is quite a generic result, independent of the details of diffusion model 
adopted for cosmic ray propagation, having its origins in such fairly tractable processes, 
such as the expected positron energy loss due to inverse Compton scattering off cosmic 
microwave background |23] and starlight. It is therefore quite unlikely that the required 



boost in DM annihilation flux can be obtained by invoking the galactic subhalo population 
— a result confirmed by detailed modelling of its effect [25]. 



While the existence of other density-borne boosts to the flux cannot, in principle, 
be completely ruled out (e.g., in addition to the subhalo population, DM caustics, tidal 
streams, mini-black holes hosting DM spikes etc. have been proposed), a large local en- 
hancement is unlikely to equally boost the positron flux at all energies, and so leave the fits 
to the flux invariant — again because of the energy dependent manner in which positrons 
lose energy as they move through the galaxy; indeed, the case for streams and caustics 



seems to have been falsified, even in principle, i.e. in terms of source function in |2f:]. The 
process of fitting the data therefore requires significant fine tuning of masses and positions 
of objects (as is also the case when invoking the standard subhalo population |^). Many 
of these proposals are, in addition, in significant confiict with, or highly constrained by, 7 
ray observations |28] |29|. 



Another class of proposals involve velocity dependent cross sections originating from 



Sommerfeld enhancement effects (e.g., |30| [32| |3^]). These can dramatically increase 
annihilation rates at velocities characteristic of galactic halos without affecting the relic 
abundance. This boost however comes at the price of fine tuning the particle masses p2|| ; 
it requires new light scalar or gauge boson, (f>, to mediate the annihilation of DM into 
leptons. The mass of this particle is constrained as: 

< am^ < few GeV, (2.6) 

where a refers to the coupling of the (p interaction squared over 47r. 

Furthermore, proposals for further enhancing the cross section for annihilation by 



considering the still lower relative particle velocities inside distant subhalos [30|, or in the 
central regions of halo cusps that are shallower than isothermal [^4| , are untenable, because, 
again, high energy positrons originating in these regions cannot safely traverse the space 
separating them from the solar neighborhood. 

In addition, Cosmic Microwave background constraints seem to rule out dramatic en- 
hancements due to low dark velocities matter velocities in subhalos, since they suggest the 
low-velocity annihilation enhancement must have saturated at the last scattering surface 
(when vdm/c ~ 10~^). WMAP results do not, on the other hand, rule out enhanced cross 
sections consistent with dark matter annihilation from the main halo that explain the cos- 
mic ray excesses. They therefore do not falsify the models described below |^5| (though 
PLANCK results may). Finally, FERMI 7-ray constraints also seem to rule out enhanced 



low particle velocity enhanced boosts in subhalos pl| [p2 |. 

3. Relic Abundance and Non-conventional Cosmology 

During the radiation epoch in the early Universe, if the DM is assumed to be in thermal 
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equilibrium, at the freeze out temperature Tp, the standard calculation for the thermal 
average of the annihilation cross section, for u <C c, yields: 



{(jv)f = a + b Vp, (3.1) 

where vp is the velocity of the DM at the freeze out temperature, which is of order 0.1 c. 
Note that the above expansion of {av) is consistent only far from s-poles and threshold. 
In our galactic halo the velocity of the DM particles is much smaller (of order 10"'^ x c). 
Its distribution is also no longer Maxwellian and its dispersion may vary significantly with 
radius. Nevertheless, due to the suppression of the kinetic energy of the DM respect to 
its rest mass, the average cross section in the galactic halo has negligible dependence on 
the velocity distribution function (we have explicitly checked this point by considering 
the isotropic distribution function associated with the NFW density profile). In this 
context, {crv)halo can be written as 

{(yv)haio ^ a. (3.2) 

This implies that the value of {crv), involved in the computation of the positron flux in 
the galactic halo, is the s-wave annihilation, which contributes to the DM relic abundance 

Given the difficulties outlined above, it seems pertinent to ask whether the constraints 
imposed on the cross section by standard relic abundance calculations are as universally 
prohibitive as they seem; and to examine mechanisms whereby they could be circumvented. 

The usual assumption is that the DM particles decoupled from the standard model 
particles when the former became non-relativistic. In that case the DM relic density is 
given by 

3l''(TF){a/TF + 3b/x],) 

where Tp is the freeze out temperature, xp = rriy-fTp ~ 20, and the number of relativistic 
degrees of freedom. It is clear that, for a ~ 10~^ GeV~'^, one gets f^^^^ ~ 10~^, which is 



inconsistent with the WMAP results |37|, as well as a host of observations concerned with 
the dynamics and large scale distribution of galaxies and clusters. 

It is nevertheless important to note that relaxing this assumption can very well give 
rise to different predictions for the relic density without affecting any of the relevant obser- 
vations. We will focus here on two such non-standard cosmological scenarios as examples 
illustrating this point, thus resolving the tension between the relic abundance constraints 
and the PAMELA and ATIC results. Our first example considers a cosmological scenario 
where the reheating temperature is associated to the decay of a standard model singlet 
field ip, in which case the DM relic density considerably increases with respect to the stan- 
dard radiation-dominated case, by virtue of the effect of direct non-thermal production by 
the ip field |^^. The second example is based on the non-conventional brane cosmology 



1 41], which also allows for several orders of magnitude increase in the cross section. The 
principal point in this first example concerns the numerical relation between the reheating 
temperature Tjih and the the DM freezing temperature Tp. The former being defined as 
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the temperature at which the oscillating field energy ceases to dominate the cosmological 
evolution, heralding the start of the radiation dominated epoch. In the standard scenario, 
where Tr// ~ 10^ GeV, the reheating epoch has no relevance in the final output of the DM 
relic density. However, for low value of Tjih, so that Trh < Tp, it can have important 



implications on the predictions of the relic abundance of DM as discussed in Ref. [38|. 

Thus, in this context, it is possible to envision the existence of an epoch in the his- 
tory of the Universe, preceding the radiation-dominated era, when the energy density was 
dominated by coherent oscillating fields; the so called reheating era. The decay width of 
this scalar field can be parameterized via 



1 m; 



2tt M2 ' 



(3.4) 



where M,, defines a high scale, acting as an effective suppression scale. 

The coupled Boltzmann equation, for the DM, scalar field tjj and radiation, is then 
solved. The detailed equations used in this computation can be found in Ref. IM]. The 



resulting relic density can then be estimated to be [39] [40| 

1 GeV^^^ /inn rio^^X 3/2 



1.5 X lO^OGeV J V {(TV 



lOOGeVV 
J 



10.75 \ 

9* J 



1/4 



(- 



VlOOGeV/ 



) ■ (3.5) 



From this equation, it can be seen that it possible to obtain ^^h?' ~ 0(0.1) with an 
annihilation cross section of order 10~^ GeV~^, as required by PAMELA and ATIC data. 
For example, for 
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Figure 3: The enhancement/suppression factor R = {VlyAi^)^/ {Vl^h'^)s as a function of the five 
dimensional scale A/5 (GeV) for = 100 (solid curve), 200 (dashed curve) and 500 GeV (dotted 
curve) . 

We now turn to the second possible scenario invoking non-standard cosmology as a 
mechanism for boosting the DM annihilation cross section. As mentioned above, it is based 
on the brain world cosmology, which is embedded in five dimensional warped space time. 
In this case, the derived Friedman equation is given by 



8ttG 



(4) 



P 1 + 



2a 



k 



+ 



(3.6) 



-7- 



where H = a/ a is the Hubble parameter and a{t) is the scale factor, p is the energy density 
of ordinary matter on the brane, while a is the brane tension; ^(4) refers to the 4Z) Newton 
coupling constant, k stands for the curvature of our three spatial dimensional and C is a 
constant of integration known as dark-radiation. This equation implies that H oc p rather 
than yjj) as in the standard cosmology. Thus, the evolution of the scale factor will be 
different from the usual one. This modification affects any relic abundance of DM that 
may be reminiscent of the radiation dominated phase of the early Universe |41] [42] [43|. 



For, in this model, the Universe undergoes a nonstandard brane cosmology at early times 
till it reaches a temperature, known as transition temperature Tt, when it sustains the 
standard cosmology. This transition temperature is defined as 

p{Tt) = 2a Tt = 0.51 x IO^^mI GeV, (3.7) 

and the transition should take place above the nucleosynthesis era (i.e., Tt>l MeV). Here 
M5 is the five dimensional Plank mass. If the freeze out temperature of the DM (Tp) is 
higher than the transition temperature, i.e., Tp > Tt and M5 < 10^ GeV, one finds that 



the ratio between the relic density in brane and standard cosmology is given by |41] 



R = {n^h\/{n^h^)s ~ 0(10^ - 10^). (3.8) 

In Fig. ^, we present the prediction for the factor i? as a function of the scale of the 
five dimensions, M5, for different values of m^, namely we consider rriy. = 100,200 and 
500 GeV. As can be seen from this figure, for M5 < 10^ the brane cosmology effect is 
quite large and the factor R becomes much larger than one. In this case the resulting relic 
density {Q.yh?)h may exceed the WMAP results Vl^h? ~ 0.1. Moreover for M5 > 5 x 10^, 
the ratio R becomes less than one and a small suppression for {Q^h?)s can be obtained. 

An analogous analysis, in the case of standard cosmology, would give {VLh?)s ~ 10~^, 
under the fixed parameters = 1 TeV , (av) = 10^^. Hence, for this calculation to be 
consistent with {Qh'^)i, ~ 0.1, R must reach at around 10^. Hence M5, in our case, has to 
be more than 10^ from the Fig. |3[ 

This brane enhancement or suppression for the dark matter relic density could be 
favored or disfavored based on the value of the relic abundance in the standard scenario. If 
($7^/i^)s is already larger than the observational limit, as in the case of bino-like particle, 
then a suppression effect would be favored and hence M5 is constrained to be larger than 
5 X 10^ GeV. However, for wino- or Higgsino-like particle where the standard computation 
usually leads to very small relic density, the enhancement effect will be favored and the 



constraint on M5 can be relaxed a bit |44]. In general, it is remarkable that in this scenario 
the dark matter relic density imposes a stringent constraint on the fundamental scale 
M5. This brane enhancement implies that stable particles with annihilation cross section 
~ 10~^ GeV~^, which seems essential for accommodating the PAMELA and ATIC results, 
remain viable DM candidates. 
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4. Extra Constraints on the Annihilation Cross Section and the Viabihty 
of our Model 



For our model to be viable, it has to be consistent with other constraints. Some of these 
were briefly mentioned earlier in our discussion; we summarize the situation as regards to 
these here. 

DM annihilation at redshift z ~ 1000 has an effect on the degrees of ionization of 
the primordial plasma, by virtue of the energy it pumps into it through the high energy 
particles it produces. By affecting its degree of ionization, it broadens the last scattering 
surface, which in turn has an effect on temperature correlations and polarization of the 
cosmic microwave background. Slatyer et. al. [35| have presented a detailed study of 
the phenomenon. Their results tend to disfavor the large boosts implied by Sommerfeld 
enhancements in dark matter substructure, but not as much the standard enhancement 
at normal local dark halo velocities. WMAP results therefore do not rule out the cross 
section boosts required for fitting the PAMELA and ATIC results in the context of non- 
standard cosmological scenarios, and are even less severe in constraining the FERMI data. 
Though further constraints from PLANCK would provide further stringent tests that may 
in principle falsify the scenarios presented here. 

Since dark matter annihilation produces secondary gamma rays, through internal 
bremsstrahlung and inverse Compton scattering, there are also severe constraints on DM 
annihilation cross sections from the observed gamma ray fluxes; in particular from recent 
observations by FERMI .The most stringent appear to be those concerning the gamma 
rays emanating dark matter dominated systems such as dwarf galaxies and galaxy clus- 
ters. Recently published results from eleven months of FERMI observations |21|, |22], seem 
to disfavor cross sections consistent with the observed cosmic ray excesses if the parti- 
cles producing these fluxes exceed the ITeV limit, particularly if the annihilation proceeds 
through the muon channel. Thus all the models presented here are (albeit in one case 
marginally) consistent with these constraints. 



Finally, the current experimental limits from CDMS [45| and XENONIO imply 
that the DM- nucleus scattering cross section is less than 10~^ GeV~^. By crossing sym- 
metry in Quantum Field Theory (QFT), the annihilation cross section and the scattering 
cross section are related; which may lead one to suspect the presence of a contradiction 
between the direct detection results and the attempt to explain the PAMELA and ATIC 
data via enhanced annihilation cross sections. But this apparent conflict can be resolved 
by assuming that the DM particle has suppressed coupling with quarks, while retaining 
reasonably large coupling with leptons; i.e., that the DM has a leptonic nature. As men- 
tioned above, this type of DM is also favored for explaining the absence of flux excess in 
the observed anti-proton flux in PAMELA, ATIC and FERMI data. It is clear that in this 
case one can naturally have the annihilation cross section into of order 10^^ GeV~^, 
whereas the scattering cross section of DM and nucleus cjdm-n is quite suppressed. 

It is also important to note that an s-wave dominated annihilation cross section of 
order 10~^ GeV~^ bears interesting corollaries concerning the nature of prospective DM 
candidates, and by implication on particle physics theory. It is well known, for exam- 
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pie, that if the DM is composed of Majorana particles, then its annihilation cross section 
into leptons is proportional to mf/m'^; which, for m-^ ~ 100 GeV, is typically less than 
10^*^ GeV~^. However, for Dirac or scalar boson type DM the annihilation cross section 
is no longer proportional to the lepton mass squared and it may thus be enhanced sig- 
nificantly. In this context, the famous DM candidate of lightest neutralino in minimal 
supersymmetric standard model, which is a Majorana particle, is not favored for explain- 
ing the PAMELA/ATIC results. On the other hand, the right-handed sneutrino, which is 
considered an interesting scalar candidate for DM in supersymmetric theories with right- 
handed neutrino can account for the large annihilation cross sections required for 
explaining the PAMELA/ATIC measurements. 

5. Conclusions 

In this paper we have proposed solutions reconciling the interpretation of the observed 
cosmic ray excesses with dark matter annihilation models by invoking non-standard cos- 
mologies. We have argued that these models have several advantages over those invoking 
density or velocity borne boost factors, especially those invoking large enhancements due 
to dark matter overdensities and halo substructures. At present, our models are consistent 
with current observations. Future data, especially those from the PLANCK concerning the 
enhanced effect of DM annihilation on the Cosmic Microwave Background, will determine 
whether our models remain viable. 
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